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We study in the framework of relativistic quantum mechanics the evolution of a system of two 
Dirac neutrinos that mix with each other and have non-vanishing magnetic moments. The dynamics 
of this system in an external magnetic field is determined by solving the Pauli-Dirac equation with 
a given initial condition. We consider first neutrino spin-flavor oscillations in a constant magnetic 
field and derive an analytical expression for the transition probability of spin- flavor conversion in the 
limit of small magnetic interactions. We then investigate ultrarelativistic neutrinos in an transver- 
sal magnetic field and derive their wave functions and transition probabilities with no limitation 
for the size of transition magnetic moments. Although we consider neutrinos, our formalism is 
straightforwardly applicable to any spin-1/2 particles. 
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I. INTRODUCTION 

In the neutrino oscillation phenomena observed so far 
oscillations occur between different neutrino flavors v\ 
(A = e, fi, t) so that an active left-handed neutrino goes 
to another active left-handed neutrino (e.g. v~ <-> v~f). 
In some situations another type of oscillations is possible, 
namely spin-oscillation, where oscillation happens be- 
tween an active left-handed neutrino and its inert right- 
handed counterpart (e.g. v\ <-> v^). Also a combination 
of these two oscillation types, so called spin- flavor oscilla- 
tions, can happen. There oscillations take place between 
an active left-handed neutrino and an inert right-handed 
neutrino of different flavor (e.g. v\ <-> v^) plQ. 

In this paper we shall consider spin-flavor oscillations 
in an external magnetic field. We assume that neu- 
trinos are Dirac particles with non-vanishing magnetic 
moments. Let us remind that Dirac particles can have 
both ordinary (diagonal) magnetic moments and transi- 
tion (non-diagonal) magnetic moments, whereas for Ma- 
jorana particles only transition magnetic moments are 
possible Q. Note that despite the recent claims of the 
experimental discovery of the Majorana nature of neu- 
trinos it is still an open question whether neutrinos 
are Dirac or Majorana particles ||. We shall use the 
formalism of relativistic quantum mechanics, that is, we 
use the Dirac equation as a starting point, which is a 
proper approach for spin-1/2 particles. One of us (MD) 
has previously used this formalism for studying neutrino 
flavor oscillations in vacuum 6] and in an external axial- 
vector field (neutrino interaction with matter). Neu- 
trino spin-flavor oscillations in electromagnetic fields of 



various configurations have been studied in Refs. [8|. The 
present paper is a continuation of these previous works. 

Let us consider a system of two Dirac neutrinos with 
non- vanishing masses and magnetic moments. In gen- 
eral, both the mass matrix and the matrix of magnetic 
moments are non-diagonal in the flavor basis of neutrino 
wave functions. When the mass part of the Hamiltonian 
is diagonalized by a unitary transformation and the origi- 
nal flavor basis is replaced by a new set of wave functions 
in the mass eigenstate basis, the matrix of magnetic mo- 
ments is generally not diagonal in the new basis. This 
means that there will be transition magnetic moments 
between the mass eigenstates. We will consider magnetic 
moment matrices of neutrinos in various bases in Sec. [ill 
In Sec. IIIII we will discuss the situation where the result- 
ing magnetic moment matrix is close to diagonal, i.e. the 
transition magnetic moments are small compared with 
the diagonal magnetic moments in the mass eigenstate 
basis, in which case one can apply the formalism devel- 
oped in Refs. [1,|3]- In Sec. lIVI we will consider the effects 
of transition magnetic moments for ultrarelativistic neu- 
trinos in a transversal magnetic field with no limitations 
on the size of any magnetic moments. In particular, we 
apply our result for studying the situation, where the 
transition magnetic moment is large compared with the 
diagonal ones. We will derive the transition probability 
for the process like vfe — > v^, where a, (3 denote two dif- 
ferent flavors in transversal magnetic field. In Sec. fVl we 
summarize our results. 



II. ELECTRODYNAMICS OF MIXED 
PARTICLE STATES POSSESSING MAGNETIC 
MOMENTS 



* Electronic address: 
tElectronic address: 



dvmaxim@cc.jyu.fi 



maalampi@cc.jyu.fi 



The magnetic moments of Dirac neutrinos are usu- 
ally non-vanishing in both flavor and mass eigenstates 
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bases [9j . This is why the flavor oscillations of neutrinos 
under the influence of an external magnetic field are as- 
sociated generally with spin-flavor conversions. One has 
to analyze the behaviour of a two Dirac neutrino system 
in a four dimensional basis 4 ,T = (V^ , v~, i/g , v^), which 
includes both chiral components of neutrinos [loj ] . The 
analytic solution of the Schrodinger evolution equation 
for this system appears to be quite complicated gener- 
ally. 

Let us denote the magnetic moments of the two Dirac 
neutrinos v ai vp as M aa , Mpp and M a p, where the last 
one is called as a transition magnetic moment. The La- 
grangian of the neutrinos in the presence of an external 
electromagnetic field = (E, B) is then given by 

C(v a , Vfj) = 22 £ (^a) _ { m !ioJ>pv a + h.c.) 

\=a,0 

~\ MxyV^VyF^. (2.1) 

AA'=a,/3 

Here, £ (^a) = v\(vy M d IJ - m\\)v\ and = 
(i/2)(7 P 7 zy — 7i/7 P ), and the parameters ttiaa and mp a 
(a vacuum mixing parameter) have dimension of mass. 

To analyze the dynamics of the system one has to set 
the initial condition by giving the initial wave functions 
of neutrinos v\ (A = a, (3) and then analytically deter- 
mine the field distributions at later moments of time. The 
analogous problem was studied in Refs. Q for mixed neu- 
trino states in vacuum and in Ref. Q for mixed neutrino 
states in an external axial-vector field. 

We assume an initial condition of the form 

i/ a (r,0) = 0, ^(r,0)=e(r), (2.2) 

where £(r) is a function to be specified. If we identify 
v a as v e or v T and up as v^, for example, this initial 
condition might correspond a situation where the source 
of neutrinos consists of pions and kaons which decay into 
muon neutrinos. 

In order to eliminate the vacuum mixing term in 
Eq. (|2.1[) we introduce a new basis of the wave functions, 
the mass eigenstate basis ip a , a = 1,2, related to the 
original flavor basis v\ through a unitary transformation 

a=l,2 

where the matrix U = (U\ a ) is parameterized in terms 
of a mixing angle 9 in the usual manner: 

/cosJJ -sin0\ 

^sinfl cos& J K ' 

The Lagrangian (|2.1[) rewritten in terms of the fields ip a 
takes the form 

C(i)x,i>2) = Y 4>Wa) 

a=l,2 

~2 £ ^a^ipbF^, (2.5) 
06=1,2 



where £o(ip a ) = ^aiyi^^n ~ maty a is the Lagrangian for 
the free neutrino ip a with the mass m a and 

Pab = Y, U aX M *X>Uy b , (2.6) 
AA'=a,/3 

is the magnetic moment matrix presented in the mass 
eigenstates basis. Using Eqs. (12.2l) - (|2.4p the initial con- 
ditions for the fermions ip a become 

V>i(r,0) = sin0£(r), ^(r,0) = cos^(r). (2.7) 

Let us assume that the magnetic field is constant, uni- 
form and directed along the z-axis, B = (0, 0,B), and 
that the electric field vanishes, E = 0. In this case we 
write down the Pauli-Dirac equation for ip a , resulting 
from Eq. (|2.5p , as follows: 

iipa=H a i> a + VTpb, a, b =1,2, a ^ b, (2.8) 

where TL a = (ap) + pzm a — p, a p^Y,^B is the Hamiltonian 
for the neutrino ip a accounting for the magnetic field, 
V = —fxpa^B describes the interaction of the transition 
magnetic moment with the external magnetic field, fi a — 
fiaa, and fi = /i a b = ^ba are elements of the matrix (/z a t,). 
Here we use the usual definitions for the Dirac matrixes 
a = 7°^, p 3 = 7° and X = 7°7 5 7. 

It is in order to remark that in Ref. the prob- 
lem of neutrino with anomalous magnetic moment was 
considered in the framework of Dirac-Schwinger equa- 
tion in an external magnetic field (Furry representation) , 
which leads to a wave equation that differs from our equa- 
tion (|2.8p . The two approaches give, however, the same 
energy spectrum [see Eq. (|A1|) ]. 

The general solution to Eq. (12. 8|) can be presented as 
follows: 

+ b«\t)vPexp(+U2Ptj\, (2.9) 

The basis spinors u a and vjf\ as well as the energy 
Ea , as a function of the particle momentum are given 
in the Appendix [X] Our main goal is to determine the 
coefficients and b a ^ consistent with both the initial 
conditions (|2.7p and the evolution equation (|2.8p . They 
are in general functions of time. 

III. PERTURB ATIVE SOLUTION 

In this section we investigate the dynamics of the 
mixed neutrinos system in the case where the transition 
magnetic moment can be considered as a small pertur- 
bation. We assume that the characteristic energy of the 
particle ip a is large compared with the energy \(iB\ associ- 
ated with the interaction induced by the transition mag- 
netic moment. Imposing the initial condition (|2.7p . we 
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now solve Eq. (|2.8[) taking the term Vipb in the Hamilto- 
nian (|2.8p as a small correction. We express the solution 
of Eq. (|2.8|) as an expansion 



^(r,t)=Vi° ) (r,t)+^ 1) (r,t) + 



(3.1) 



where ipi^ (r, i) is the solution of Eq. (|2.8[) when V^ft = 
0, that is, the eigenvalue of the Hamiltonian H a - The 
function ip$p (r, t) is linear in /j,B, and the ellipses stands 
for terms of higher order in fiB. 

Let us first search for the zeroth order solution ip a v 
The coefficients a„ and &i defined in Eq. (|2.9[) are for 
^4 time independent. Using the orthonormality con- 
ditions for the spinors and Va [see Eq. (|A4I) ] and 
following the method developed in Refs. @, we find 



do) 



^°)(r,i) 
where 

^a(p,<) = 



e ipr S Q (p,^a(p,0), (3.2) 



d 3 p 

(2tt) 3 

£ [(«<«®u«)t)«p(-^«) t) 
±i 



C=±i 

+ 



exp (+i^ c) t) 



(3.3) 



is an analog of the Pauli- Jordan function for a spin-1/2 
particle in an external magnetic field and 



d 3 re- ipr V> a (r, 0) 



is the Fourier transform of the initial wave function (12. 7| . 

Eqs. ((S2D and together with Eqs. (13]) and (j2TT|) 
[or Eq. (|2.2p ] allow one to describe the evolution of any 
system of two Dirac fermions in an external magnetic 
field. These expressions directly follow from the Lorentz 
invariant Pauli-Dirac equation (|2.8|) . and they are valid 
for an arbitrary initial condition. In particular, the ini- 
tial momentum of the particle can be arbitrary, though 
neutrinos are of course generally ultrarelativistic. 

To proceed we have to specify our initial condition, 
that is, give the function £(r). Following our previous 
studies [6J, we choose the initial condition for vp as a 
plane wave £ (r) = e lkr £o , with the initial momentum be- 
ing directed along the :r-axis, i.e. k = (fc, 0, 0). Given 
that we have chosen the magnetic field as B = (0, 0, B), 
we are thus studying the evolution of neutrino wavefunc- 
tion in a transverse magnetic field. We also consider neu- 
trinos as ultrarelativistic particles, i.e. k 3> toi,2, imply- 
ing £q = (1/2)(1, — 1, — 1, 1). It is easy to see that £o 
obeys the condition (1/2)(1 — £i)£o=£o and is normal- 
ized to one. Hence, the spinor £(r) describes an ultra- 
relativistic particle propagating along the x-axis with its 
spin directed opposite to the x-axis, i.e. a left-handedly 
polarized neutrino. 

In contrast to the mixing due to a mass matrix, which 
is always present, the magnetic transition only occurs in 
the presence of a magnetic field. Let us assume that a 



particle passes through a region of length ro where there 
is a nonvanishing magnetic field, i.e. B^0if0<t<to 
and B = if t < and t > to where to ~ r$. The 
particle emission and detection are assumed to take place 
while the magnetic field is switched off, at t e = — St 
and td = to + St, respectively, where St is some small 
period of time. Hence, the operator Ei describes the 
polarization state of the initial and final particle, and 
one can assume that the system is initially prepared to a 
state described in terms of the spinor £o as we have done 
above. Of course, the wave function must be continuous 
in the borderlines of the B — and B ^ regions. 

One can now determine the field distribution of the 
particle v a , which was initially absent. With help of 
Eqs. (USD, (CUD and O, and by taking into ac- 

count the initial condition we have chosen, we obtain the 
following expression for the wavefunction of the right- 
polarized state of v a : 



v£ )R (x,t) =£(l + £i) cos 0^1"' (s,t)- sin fl^fot) 



1 



(0), 



',(0)/ 



2 

=i sin 9 cos 9 [e~ i£lt sm((j,iBt) 



AS 2 t 



sin(/x 2 Bt)\e ikx K , 



(3.4) 



where 



+ k 2 , 



T 



(1/2)(1, 1,1,1). (3.5) 



We have used here the identities 



(« ± ®« ± t)f = 0, 
i(l + Ei)(u ± «u ± t)&=±^K 0j 



(3.6) 
(3.7) 



which result from Eq. (|A6[) . Note that in the ultrarel- 
ativistic limit m a /k <C 1 one can neglect the neutrino 
mass dependence of the basis spinors and and 
this is why the subscripts a and b are omitted in Eq. (13. 7p 
(see also Sec. IIVI and Eq. (|A6j) for the explicit forms of 
u^' and v^'). It should be mentioned that Eq. (|3 . 6|) 
implies that no particles with negative energies (antipar- 
ticles) are produced in neutrino interacting with an ex- 
ternal magnetic field. 

The measurable quantity is the square of the wave 
function Va (x, t), which can be interpreted as the tran- 
sition probability of the process v\ — > v^. From Eq. 
we obtain 



where 



,(0)R 



sin (29) { sin 2 (SiiBt) cos 2 (pBt) 
sin(Mi-Bt) sm(fi 2 Bt) sin 2 [$(fc)i] }, (3.8) 



$(fc) 



Sm 2 



(3.9) 



is the phase of vacuum oscillations, Sin 2 
5fX = (fit - u 2 )/2 and p,= ((i 1 + n 2 )/2. 
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In the case of neutrinos having equal magnetic mo- 
ments, Hi — fi2 = Mo 7 Eq. (|3.8[) leads to the result 
one would expect. Namely, Eq. (|3.8p can be rewrit- 
ten as P = P F Ps, where P F = sin 2 (26 l ) sm 2 {$(k)t] is 
the usual transition probability of flavor oscillation and 
Ps = sin 2 (/io£i) is the probability of the transition 
between different polarization states within each mass 
eigenstate. That is, as the magnetic moment interac- 
tions are in this case insensitive to flavor, the transitions 
between flavors are solely due to the mass mixing. 

Let us now consider the first order correction tp^](r, t) 
in Eq. (|3.ip . Using the same method as in Ref. [3] we 
obtain 



V>i 1} M) = -i 



d 3 P 



3 ipr 



C=±i 



(2tt) 3 

x exp {-iE^H)Vg { P + (vP ® w<«t 
x exp(-Hi#>t)U7^Wp,0), (3.10) 



where 



gM= / dt'e X p(+iEp t ')S b (p,t'), 
Jo 

Jo 



In Eqs. (gUDl) and (j3~TTj) a ^ 6. 

With help of Eqs. ([2"3]) . (j^4|) . (f3TT0|) and (|3~TT|) the 

first-order correction to z'i '' gets the form 



v^ K {x,t) =i cos 2(9 



, iSt sinAt iCTt sin<5i ^ Jfca . 

^ I " 1 i~ r 



(3-12) 



where 



S = £ + p,B, A = 5£ + SfiB, 

a = £-pB, 5 = 5£-6fiB, (3.13) 



and 



66 = 



Si- £2 



£ = 



£i+£-2 



(3.14) 



It can be seen from Eqs. (|3.I2p and (|3.13|) that the addi- 
tional constraint should be fulfilled for the perturbative 
approach to be valid, namely S£ ^ ±5{iB. It should be 
noticed that Eq. (|3.I2|) is obtained in the ultrarelativistic 
limit, where small terms m a /k <C 1 are neglected. 

Let us now find out the effect of the first-order cor- 
rection to the transition probability. Using Eqs. 



and (|3.12p we obtain 

P^n(t) =4 0)Rt -i 1)R + h.c 



(3.15) 



sin AO 



4[$ 2 (fc) - (6fj,B) 2 ] 
x j$(fc)sin[2$(fc)t] sin(2c^Bi) 
-4(5 fiB ( sin 2 (6fj,Bt) cos 2 {jlBt) 
+ sm(fiiBt) sin(^ 2 Si) sin 2 ) }. 

The transition probability of spin-flavor oscillations be- 
tween mass eigenstate neutrinos up to effects linear in 
the transition magnetic moment is given as a sum of the 
probabilities given in Eqs. (|3.8j> and (|3. 15[) . 

It should be noticed that the method used above allows 
one to examine particles with arbitrary initial distribu- 
tions £(r) since Eqs. and (O, as well as Eqs. (f3~TU|) 
and (|3.Iip . are not restricted to any specific initial con- 
dition. Also, as we mentioned before, the initial mo- 
menta of particles is not restricted. This has relevance, 
of course, only in the case the results are applied to other 
particles than neutrinos. Nevertheless, our analysis is not 
totally general as we have assumed the transition mag- 
netic moment small. It turns out, actually, that one can 
circumvent this restriction in the case of ultrarelativis- 
tic particles. That is, in the case of neutrinos one can 
solve, as we will do in the next Section, the Pauli-Dirac 
equation analytically for an arbitrary magnetic moment 
matrix. 



IV. EVOLUTION OF ULTRARELATIVISTIC 
PARTICLES WITH ARBITRARY MAGNETIC 
MOMENTS MATRIX 

In this section we study the influence of the transition 
magnetic /x on the evolution of the mixed neutrino system 
without assuming fi to be small. 

Using the orthonormality conditions of the basis 
spinors [see Eq. (j A4|) ] and the fact that ±Ei^ are the 
eigenvalues of the Hamiltonian H a , a = 1,2, 

U a u{0 = +40.(0, HaV (0 = - E p v (O t (4.!) 

one obtains from Eq. (|2.9p the following ordinary differ- 
ential equations for a^\t) and ba (t): 

iflW =e X p(+iEPt)upV 

x J2 [aPuPe Xp (-i E Pt) 
C'=±i 

+ bpvPe*p(+iEpt)] ) 
lb® =exp(-LE(^H C)t ^ 

x E [apupexp(-iEpt) 



C'=±i 



bl C) vl c> exp(+iEpt)]. 



(4.2) 
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Eqs. (|4.2|) are subject to the initial conditions 
^ C)(0)= (2^ U P^ a(P ' 0) ' 

^ c) (o)=^^(p,o), 



(4.3) 



which result from Eq. (|2.9[) . Note that the functions a& 
and are here time dependent in general, in contrast 

to sec. inn 

Let us choose the initial wave function of vp [see 
Eq. (|2.2p ] as £(r) = e lkr £o, with the initial momentum 
being aligned along the x-axis, k = (k, 0, 0). The mag- 
netic field is assumed to be B = (0, 0, B), implying that 
we study the propagation of neutrinos in the transver- 
sal magnetic field. Given that k _L B, we can rewrite 
the Eq. (|4.2[) for the function ai in the form (the corre- 
sponding equation for the function can be obtained 
analogously) 



at(ut\V\ut)exp [i(£± - E±)t] 
+ tf«\VK)exp[i(Et + E*)t}. 



(4.4) 



Here we use the explicit form of the basis spinors in the 
transversal magnetic field given in Eq. (|A2[) . Note that 
(ut\V\uf) = and (ut\V\v±)=0. 

A further simplification of Eqs. (|4.2p is obtained when 
we study the ultrarelativistic initial wave function, k 3> 
mi, 2l and = (1/2)(1, -1, -1, 1) (see Sec. |m|, in 
other words the system is, like in the case we studied 
in the previous Section, in the state initially. With 

help of the obvious identities (u^;\V\uf) = T^B and 
(m^IV^u^) = 0, which result from Eq. (|A6|) . one can cast 
Eq. (|4.4|) into the form 



= Taf f iBexp [i(Et-E±)t 



(4.5) 



Let us note that the Eq. (|4.5p is similar to the evolution 
equations for a neutrino interacting with a twisting mag- 
netic field, which were examined in Refs. On the 
basis of this previous study we are able to write down 
the solutions as 



af(t) =F ± af{0) + G ± af(0), 
af(t) =F ± *af(Q)-G ± *af(Q), 



(4.6) 



where 



F ± 
G ± 



cos Q±t — i— — — sin il±t 



2Vt A 



exp (iw±t/2), 



± i— — sin Q±t exp (icj±t/2), 



(4.7) 



and 



0± = y/ifiB) 2 + (^±/2) 2 , uj± = Ef - Ef. (4.8) 

The derivation of Eqs. (|4.6|) - (|4.8|) from Eqs. (|4.5[) is pre- 
sented in Appendix 151 



Using Eq. (|2T9|) and Eqs. (|42| - (|4~8| and the identity 
<g> i)K)t) £ = [see Eq. (|377|) ] we obtain the wave 
functions ip a , a = 1, 2, as, 

ipi(x, t) = exp (-i-E^t) (u + (g) u +t ) 

x [F+VxMHG+^M)] 
+ exp (-iEft) (u~ <g> u~ t ) 
x [F-ih{x,0) + G-fo{x,0)], 

ijj 2 (x, t) = exp (-iEft) (u + <x> u +t ) 

x [F + *V 2 (a; ! 0)-G + *Vi(a;,0)] 

+ exp (— iE^t) (w~ ® u -t ) 

x [F-*ifo{x, 0) - CT Vifo 0)], (4.9) 

which satisfy the chosen initial condition since G ± (0) = 
0, F ± (0) = 1 [see Eqs. ([47T]) ] and [(u+ g> u+t) + (u" ® 
u-^)]ip a (x,0) = 4> a {x,0) [see Eq. ([A"6])]. Note that the 
subscripts a and 6 are again omitted in the basis spinors 
and as we assume ultrarelativistic particles. 
With help of Eqs. (US}, (HP and gl]) [see also 
Eq. (|3.4p ] we receive for the right-handedly polarized 
component of v a the expression 



v ^{x, t) =- { sin 6 cos 6 [e~' l£lt { e 1 ^ Bt F+ 



e -^ lBt F-) 

e ~ i£ 2t ^ViBt p+* _ e -iM2-Bt^-*^j 

cos 2 9e' i£lt (e if " lBt G + - e-^ lBt G~) 



sin 2 6e~ i£2t {e it22Bt G + * - e -^ 2Bt G-*)} 



x e lkx n . 



(4.10) 



The quantities £ a and kq are defined in Eq. (|3.5p . We 
have used here the identity given in Eq. (|3.7p . 

Finally, taking into account Eqs. (|4.7p and (|4.8p it is 
possible to express the wave function in Eq. (|4.10p in the 
form 



v*{x,t) 



sin 9 cos 6* — 
2i 



— — sin(f2+i) exp (i/iBt) 

1 



^— sin(f2_i) exp (—i/iBt) 



■LxB 



x cos(fiBt) > exp (— i£i + ifco;)^- 



sin(f2_t) o „ 
■ sin 9 



(4.11) 



can 



The transition probability for the process 
be directly obtained as the squared modulus of v^{x,t) 
from Eq. (|4~10|) or Eq. (|4~TTj) . that is P„l_„,r(*) = 

\v^(x,t)\ 2 . Notice that the probability is a function of 
time alone with no dependence on spatial coordinates. 
This is of course obvious as we have taken the initial 
wave function as a plane wave and the the magnetic field 
spatially constant. 

Let us now apply the general results Eq. (|4.10p or 
Eq. (|4.1ip to two special cases. We first consider the sit- 
uation where nx,2 3> i-e. the case when the transition 
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magnetic moment is small compared with the diagonal 
ones. Using Eqs. (|4.7[) and (|4.8p we find that in this case 
F ± w 1 and tt± w w±/2, and Eq. ([4~T0]) takes the form 



i< smf cosf 



\e~ i£lt sin nxBt - e 



-iS 2 t 



sin /i2-B^l 



cos 26> 



^5 



sin Ai 



sin St 



x e ifex K - 



(4.12) 



In the previous Section we studied this case perturba- 
tively, and one can easily check that results obtained 
there [see Eqs. (pRj) and (|3~T2"]) ] coincide with P~T2"]) . 

As another application of our general result we will 
study the situation, where the transition magnetic mo- 
ments are much larger than the diagonal ones, that is 
M ^ Mi, 2- In this case Eqs. (|4.7[) gives F + « F~ and 
G+ w — G~, and we receive from Eq. (|4. 10[) or Eq. (|4.11|) 
for the wave function the expression 



v%{x, t) ~i exp {—i£t + ifcx) 

x cos(2#)— sin(Ot)/«o, 



where 



n = V(^) 2 + $ 2 (fc). 



The transition probability for the process i/g 
then given by 



Pu^{t) - cos 2 (20) 



^5 

7T 



sin 2 (fti). 



(4.13) 



(4.14) 



is 



(4.15) 



The behavior of the system in this case is schematically 
illustrated in Fig. [T] It should be noticed that the analog 
of Eq. (|4.15[) was obtained in Ref. [|J where the authors 
studied the resonant spin-flavor precession of Dirac and 
Majorana neutrinos in matter under the influence of an 
external magnetic field in frames of the quantum mechan- 
ical approach. 

Let us discuss the applicability of our results to one 

specific oscillation channel, ^ v^. According to the 
recent experimental data (see, e.g., Ref. [l2|) the mixing 
angle between and v T is close to its maximal value 
of 7r/4. In the limit of maximal mixing the magnetic 
moment matrix given in Eq. (|2.6|) takes the form 



(Mob) ~ 

(M TT + M w )/2 + M Til -(M TT - 

-{M TT - M w )/2 (M TT + M w )/2 - M Tii 



(4.16) 



Eqs. d3HI) and (j3~T5j) [or Eq. gH])] are valid in 
the case this matrix is close to diagonal, i.e. when 
\{M TT - M w )/2| « \{M TT + M w )/2 ± M Tll \. In con- 
trast to the mixing angles, the experimental information 
about the neutrino magnetic moments is very limited. 




FIG. 1: The schematic illustration of the system evolution in 
the case [i S> fli,2- The horizontal lines of the figure corre- 
spond to various neutrino eigenstates at different moments of 
time (t = and i). The expressions next to arrows corre- 
spond to the appropriate factors in the formula (|4.13fl of the 
wave function. The arrows from Vg{Q) to (0) and ip2 (0), 
for example, indicate the vacuum mixing matrix transforma- 
tion at t = and the arrow from ipi (0) to ip2(t) the evo- 
lution of the mass eigenstates with the helicity change. The 
transitions ^a(O) — > ipt'it) can be described by the formula, 
ftW = Al^r^^O), where A L ^ n = i{i±B/Q) smQ,t. 



Also theoretically, very little is known about their mag- 
nitude. What one knows is that the diagonal magnetic 
moments M\\ could be very small in the extensions of 
the Standard Model M\\ ~ f CP 19 (mAA/eV)/iB (see, e.g., 
Refs. [Hj|). The contributions to the transition magnetic 
moments, M TM in our case, can be much bigger close to 
the experimental upper limit of 10 _10 /iB [14] . One can 
see that for any conceivable values of the masses of the 
known neutrinos, and M TT are orders of magnitude 
smaller than 10 _10 /xb- Our results [Eq. (|3.8[) ] are appli- 
cable in this kind of situation. 

We can also consider the opposite situation of M\\ ^> 
M Tfl . It is implemented, e. g., i n the minimally extended 
standard model (see Refs. [l3j) when the magnetic mo- 
ments matrix is diagonal in the flavor basis, i.e. M Tfl is 
negligible. If it is the case, one observes that [i\ — fj,2 in 
Eq. (|4.16p . Thus uj + = lu- = lu and £1 + = fi_ = fl [see 
Eqs. flUD and fA"5j) ]. Eq. ([4~TT|) takes the form, 



<{x,t) = —sin(Qt)sin{fiBt) 



2fl 

x exp ( 



-if t + ikx)no, 



(4.17) 



where we account for that cos 9 = sin 6 = 1/V2. The 
transition probability can be calculated on the basis of 
Eq. (HTTP . 



$(fc) 
n 



sin (fit) sin 2 (fiBt). (4.18) 



Now the magnetic moments are (x = (M TT + M MAt )/2 
and \i = -(M TT - M w )/2. In Eq. ([435]) the phase of 
vacuum oscillations, $(fc), is determined in Eq. (|3.9[) and 
the parameter ft is introduced in Eq. (|4.14|) . 

At the end of this section it should be noticed that 
the our results for the description of spin-flavor oscil- 
lations of relativistic Dirac neutrinos are consistent with 
the standard quantum mechanical approach based on the 
Schrodinger equation of evolution. The consistency is 
briefly discussed in Appendix [Cl 
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V. SUMMARY 

We have studied in the framework of relativistic quan- 
tum mechanics the evolution of a system of two Dirac 
neutrinos that mix with each other and have non- 
vanishing magnetic moments. By solving the Pauli-Dirac 
equation with a given initial condition we determined the 
time evolution of this system in an external magnetic 
field. We applied the recently developed approach for 
the description of the first quantized mixed particles (see 
Refs. @,[|). 

We first (Sec. IIIip studied the special case of the mag- 
netic moments matrix, which is close to diagonal in the 
mass eigenstates basis. We obtained the general expres- 
sions for the fermion fields distributions exactly account- 
ing for an external magnetic field in the leading and the 
next-to-leading orders in the transition magnetic moment 
in the mass eigenstates basis. As we started with the 
Lorentz invariant neutrino wave equations, the derived 
neutrino wave functions are valid for arbitrary initial con- 
ditions. For instance, with help of Eqs. (|3.2[) and (|3.10[) 
one can describe the propagation of non-relativistic par- 
ticles in an external magnetic field. 

We considered spin-flavor oscillations of ultrarelativis- 
tic neutrinos in a constant transversal magnetic field and 
derived an analytical expression for the transition proba- 
bility for the spin-flavor conversion in the limit where the 
energy associated with the magnetic interaction arising 
from transition magnetic moment is small compared with 
the total neutrino energy. Then (Sec. IIV[) we examined 
the evolution of ultrarelativistic Dirac neutrinos in the 
constant transversal magnetic field. In this situation it 
was possible to take into account the influence of all mag- 
netic moments on the neutrinos evolution, in contrast to 
the analysis in Sec. IIIII where the effect of transition mag- 
netic moments where considered as a small perturbation. 

Our results are valid for arbitrary magnetic field 
strength. The general analytical formula (|4.11[) for the 
neutrino wave function allows one to calculate the tran- 
sition probability and consider various limiting cases. To 
our knowledge, the general case of spin-flavor oscillations 
of Dirac neutrinos with an arbitrary magnetic moment 
matrix has not been studied analytically before. Both 
experimental and theoretical information about magnetic 
moments of Dirac neutrinos is still very limited (see, i.e. 
Ref. [Hj]). In this light our results may turn out useful as 
they allow one to describe phenomenologically spin-flavor 
oscillations of Dirac neutrinos in any model which pre- 
dicts neutrino magnetic moments. Although we consider 
neutrinos, our formalism is straightforwardly applicable 
to any spin-1/2 particles. 
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APPENDIX A: SOLUTION TO THE 
PAULI-DIRAC EQUATION 

The general solution to the Pauli-Dirac equation of the 
form [see also Eq. (|2.8|) ]. iip a = Haipa, for a massive neu- 
tral fermion with the anomalous magnetic moment was 
found in Ref. [l5[ ■ The energy as a function of the par- 
ticle momentum p = (pi,P2,P3) and the spin quantum 
number £ = ±1, which characterizes the spin direction 
with respect to the magnetic field, has the form 



£P=£ a -Cn a B, £ a = sjml+p 2 1 +pl (Al) 
The basis spinors are expressed in following way: 



,(0 



(0 \ 4>ioi a 



1 



2VK 



(o i -<f>i<xi 



(A2) 



(A3) 



where 



(t>t = y/1 ± (m a /£a, at = yE ( a Q ± , 

and tamp — Pz/pi- The quantum number £ = ±1 char- 
acterizes the direction of the particle spin with respect 
to the magnetic field since the polarization operator 

II Q = m a S 3 + p 2 [S x p] 3 - p, a B, 

where p 2 = i7°7 5 , commutes with the Hamiltonian 7i a . 
We also present the orthonormality conditions for the 
spinors and vi^, 

uPup = v^vp = 5 CC , uP v p = 0, (A4) 

which follow from Eq. (|A2j) . The details of the derivation 
of Eqs. dXT|) and (J5J) can be found in Ref. 

If we study the propagation of ultrarelativistic particles 
in a transversal magnetic field, Eqs. (|A1[) and (|A2[) take 
the form 



=£ a - Qi a B « p + -2. - Cti a B, (A5) 
2p 
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and 



,+ -. 



1 

: 71 







'0\ 





1 


1 




u — —= 
V2 


1 


















I] 




(i 









(A6) 



In Eqs. (|A5|) and (|A6|) we assume that the momentum 
is directed along the x-axis, p = (p, 0, 0). The subscript 
a is omitted in Eq. (IA6j) since we neglect small terms 
(rn a /p) <C 1 there. 



APPENDIX B: SOLUTION TO THE ORDINARY 
DIFFERENTIAL EQUATIONS FOR THE 



FUNCTIONS ai° AND b% 

Let us study the time evolution of the two compo- 
nent spinor Z T = [Z\,Z<z) which is governed by the 
Schrodinger equation of the form, 

iZ = HZ, (Bl) 

where the Hamiltonian has the following form: 

H = (MBiy iwt e Q . (B2) 



Here fj,, B and lo are the real parameters. Eq. (|B1[) should 
be supplied with the initial condition Z(0). To find the 
solution to Eqs. (|B1[) and (|B2|) we introduce the new 
spinor Z' by the relation, Z = WL 1 ', where the unitary 
matrix hi reads 



U 



ot/2 







e' 1 ^ 2 



Now Eq. (|B1[) is rewritten in the following way: 
iZ' = H' Z', 



(B3) 



(B4) 



with the new Hamiltonian H ' which is obtained with help 
of Eqs. (TBTll-(rB3ll. 



n> = l*HU-*M=^g _f /2 ) (B5) 

Note that the initial condition for the spinor Z'(0) is the 
same as for Z(0), Z'(0) = Z(0), due to the special form 
of the matrix hi in Eq. (|B3[) . 

Supposing that the Hamiltonian H' in Eqs. (IB4|) 
and (|B5[) does not depend on time we get the solution 
to Eq. dEl) as 

Z'(t) = exp(-Lff'i)Z / (0) 

= (cos/2t-i(crn)sinJ7t)Z'(0), (B6) 



where n = (hB,0,lo/2)/(1 is the unit vector, ft = 
\J {\iB) 2 + {oj/2) 2 and a are the Pauli matrices. Using 
Eqs. (|B3[) and (|B6|) we arrive to the expressions for the 
components of Z written in terms of the initial condition 
Z(0): 



= (cos 



i^sin(^)e i " t / 2 Z 2 (0), 



Z 2 (t) = 



cos fit 

77 



l- 



t/2 Z 2 (0) 

sin(^)e- i " t / 2 Z 1 (0). (B7) 



Then we identify the spinor Z with different functions 
and fri'' •* , the parameters B and lo - with the strength 
of the magnetic field and energy differences respectively 
(of course, with the proper signs). With help of Eqs. (|4. 5f> 

(and analogous equations for the functions b^) one gets 
the list of four cases: 

• For Z T = (a±, a^), B = —B, lo = uj + and Q = Q + ; 

• For Z T = (a^, a^), B — B, lo — u>- and Q = f2_; 

• For Z T = B = B, lo = -lo + and fl = 

• For Z T = (b^,b-), B = -B, lo = and fl = 

Using these formulae together with Eqs. (|B7[) we readily 
arrive to Eqs. (|4.6[) - (|4.8[) . Note that the dynamics of 
the system (JbTJ) and (|B2() is analogous to the quantum 
mechanical description of neutrino spin-flavor oscillation 
in a twisting magnetic field studied in Refs. pTj |. 



APPENDIX C: SCHRODINGER DESCRIPTION 
OF SPIN-FLAVOR OSCILLATIONS 

Let us study the Schrodinger evolution equation for the 
two mass eigenstates neutrinos with magnetic moments 
in an external transversal magnetic field, 

i* 



H = 




(CI) 



where mass eigenstates energies £12 are given in 
Eq. (|3.5p . The wave function of neutrinos should be pre- 
sented in the following form 10J: ^f? T = (4>i,tp%, ipi, V^); 
where V'i 2 are one-component objects. The initial con- 
dition consistent with Eq. (|2.7p is 



(C2) 



f T (0) = (sin 6, cos 9, 0,0). 



To study the quantum mechanical evolution of the 
system we look for a solution of Eq. (|C1|) of the form 
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~ e lAt . Solving the corresponding secular equation 
we get the eigenvalues of the Hamiltonian H in the form, 

A+ 2 = £ + ±ft+, ai; 2 = s_±o_, (C3) 

where E± = [E^+E%)/2 and fi± are given in Eq. (I4.8|) . 
The energy levels in an external magnetic field £^ are 
presented in Eq. (jA5|) . The general solution to Eq. (jClj) 
is thus takes the form 



#(t) =(aiwie- i0 +* + a2U2e ln +*)e- ls+t 
+ (/3iWie- i0 -*+/3 2 w 2 e 

where the basis spinors are given as 



(C4) 



Ml, 2 



«1,2 



\t»b/r± 

(TVB/S T \ 

V 5 T 



2,/fT 



(C5) 



where R± = v / n+±u!+/2 and S± = y/fl- ±u)-/2. The 
spinors Uip and V1.2 in Eq. (|C5[) are the eigenvectors of 



the Hamiltonian 7f corresponding to the eigenvalues A^ 2 

and A^~ 2 i respectively. 

The coefficients ai,2 and Bip should be chosen so that 
the initial condition in Eq. (|C2[) is satisfied. We choose 
them as 



Oh 



2R ± 



sin 6 1 ± 



2^4 



T 7— cos e 

i fi-L 



(C6) 



Bx 2 = i sm6 1 ± ± — cos( 



Then, using Eq. (|C4[) . we arrive to the following helicity 
components of the wave function 'J: 

^(t) = I|( e - i£ i +i F+ + e- iE i t F-)sm6 
+ (e- iE i t G+ +e- i ^*G-)cos^}, 

V^(<) = I|( e - i£ 2 + tF+* + e - iE i t F-*)cos6 
_ ( e -i£+t G +* +6-^*0-*) sin 6»|, 

(t) =1 1 ( e - i£ i +t F + - e- iE ^ t F-)sm 9 
+ (e- iE i t G+ - e - i ^*G-)cos6i}, 



2 R w 4{( e_1 



S 2 t F +*_ e -iB 2 ^ F -*) cos( 



e 



k )sin6i}, (C7) 



where the functions f ^ and G are introduced in 
Eq. (@2J). 

Let us compare these wave functions with the results of 
our approach, presented in the main text. Initial condi- 
tions for the four-component mass eigenstates wave func- 
tions are tpi(x,0) = sin#e lfca: £o and ip2[x, 0) = cos#e lfcE £o 
[see Eq. (12. 7p ], where the spinor £0 was introduced in 
Sec. IIIII It is easy to check that 



^i,a(a. °) =2 f 1 - s i)^i.2(^, 0) = ^(z, 0), 



<aM) =i(i + 20^1,2(3;, 0) = 0. 



These expressions imply that no right-polarized neutrinos 
exist initially. 

With help of the following identity: 



^(l-EiHt^Sti^Ho^o, 



which can be derived from Eq. (|A6I) . as well as using 
Eqs. (|3.7[) and (|4.9[) , we obtain the helicity eigenstates of 
the four-component wave functions ip a , a = 1,2, 



+ 1 
+ ' 



e -^i + ti?+ + e -iS 1 t J? -) sin6) ( Cg j 
; -i-B 1 + *G+ + e -iS 1 -t G -) cos ^J e ife^ 0) 



-iE. 



*G+* + e^" ( G~*) sm6>}e i,ca: £ 
-i£rt F -) sin6 , 

i£ i +t G+ - e- i ^ t G-)cos0}e ifcx K O , 



e - 1B + tF+ * _ e -iE 2 t F -*j CQs0 



«0i 



where the spinor «o is given in Eq. (|3.5j) . The compar- 
ison of Eq. (|C8j) with Eq. rtC7[) shows that Schrodinger 
approach gives results analogous to the results obtained 
in our approach based on Pauli-Dirac equation. 



10 



[1] J. Schechter and J. W. F. Valle, Phys. Rev. D 24, 1883 
(1981); M. B. Voloshin, M. I. Vysotskii, and L. B. Okun', 
JETP 64, 446 (1986); Sov. J. Nucl. Phys. 44, 440 (1986); 
E. Akhmedov, Phys. Lett. B 213, 64 (1988). 

[2] C.-S. Lim and W. J. Marciano, Phys. Rev. D 37, 1368 
(1988). 

[3] See, e.g., M. Fukugita and T. Yanagida, Physics of neu- 
trinos and applications to astrophysics, Springer, Berlin, 
2003, p. 593. 

[4] H. V. Klapdor-Kleingrothaus, et al, Phys. Lett. B 586, 
198 (2004), hep-ph/0404088. 

[5] S. R. Elliott and J. Engel, J. Phys. G 30, R183 (2004), 
hep-ph/0405078. 

[6] M. Dvornikov, Phys. Lett. B 610, 262 (2005), hep- 
ph/0411101; in Proceedings of the IPM school and con- 
ference on Lepton and Hadron Physics, Tehran, 2006, ed. 
by Y. Farzan, eConf C0605151 (2007), hep-ph/0609139; 
hep-ph/0610047. 

[7] M. Dvornikov, Eur. Phys. J. C 47, 437 (2006), hep- 
ph/0601156. 

[8] A. Egorov, A. Lobanov, and A. Studenikin, Phys. Lett. 
B 491, 137 (2000), hep-ph/9910476; M. S. Dvornikov 



and A. I. Studenikin, Phys. At. Nucl. 64, 1624 (2001); 
ibid. 67, 719 (2004); M. Dvornikov, ibid. 70, 342 (2007), 
hep-ph/0410152. 
[9] K. S. Babu and V. S. Mathur, Phys. Lett. B 196, 218 
(1987). 

[10] W. Grimus and T. Scharnagl, Mod. Phys. Lett. A 8, 1943 
(1993). 

[11] A. V. Borisov, V. Ch. Zhukovsky, and A. I. Ternov, Sov. 

Phys. J. 31, 228 (1988). 
[12] M. Maltoni, T. Schwetz, M. A. Tortola, and J. W. F. 

Valle, New J. Phys. 6, 122 (2004), hep-ph/0405172. 
[13] K. Fujikawa and R. E. Shrock, Phys. Rev. Lett. 45, 963 

(1980); M. Dvornikov and A. Studenikin, Phys. Rev. D 

69, 073001 (2004), hep-ph/0305206. 
[14] W.-M. Yao et al., J. Phys. G 33, 1 (2006). 
[15] I. M. Ternov, V. G. Bagrov, and A. M. Khapaev, JETP 

21, 613 (1965). 

[16] A. Yu. Smirnov, Phys. Lett B 260, 161 (1991); E. Kh. 
Akhmedov, P. I. Krastev, and A. Yu. Smirnov, Z. Phys. 
C 48, 701 (1991). 



